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(Dated: June 15, 2011) 

Abstract 

We present the results of a Dalitz plot analysis of D° — > Kg7T°7T using the CLEO-c data set 
of 818 pb _1 of e + e~ collisions accumulated at y/s = 3.77 GeV. This corresponds to three million 
pairs from which we select 1,259 tagged candidates with a background of 7.5 ± 0.9 percent. 
Several models have been explored, all of which include the K*(892), K 2 *(1430), if* (1680), the 
/o(980), and the cr(500). We find that the combined 7r°7r° S-wave contribution to our preferred fit 
is (28.9±6.3±3.1)% of the total decay rate while D° -> K*(892)°vr contributes (65.6±5.3±2.5)%. 
Using three tag modes and correcting for quantum correlations we measure the D° — > Kg7r°7T 
branching fraction to be (1.059 ± 0.038 ± 0.061)%. 



*Now at: Pacific Northwest National Laboratory, Richland, WA 99352 
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I. INTRODUCTION 



The substructure of D° — > K^h + ti~ has been the object of intense recent study due to 
its relevance to CKM physics but there is currently relatively little information on the 
substructure of D° — > Kg7c°7r°. A 1993 CLEO II publication based on 206 events found the 
branching fraction of the decay chain D° — > K*(892)(— > K^ir )^ to be (6.7±i:|l x 10" 3 and 
the branching fraction of the non-resonant contribution to be (4.5 ±1.1) x 10 -3 [2j. The total 
D° — > KgTT°7c° branching fraction has recently been measured by the CLEO collaboration 
to be (8.34 ± 0.45 ± 0.42) x 10 -3 in an analysis that used 281 pb" 1 of data at = 3.77 
GeV 0|. 

The pursuit of a more comprehensive study of the substructure of D° — > K^^tt is 
motivated in part by the fact that this mode should provide a cleaner way to investigate 
the 7T7T S-wave substructure observed in D° — > Kg7i + 7T^ 0-0] ■ Analyses performed by both 
BaBar and Belle use eight nir resonances, four of which are spin- zero (including two distinct 
a resonances , to fit the D° ->• K%it + it- Dalitz plot, but the -kit S-wave components in 
these analyses are masked by the very sizeable P-wave contributions of the p resonance. 
By studying K^tt°tt we eliminate final states involving p mesons and expect any S-wave 
structure present in the decay to be more prominent. 



II. EXPERIMENTAL TECHNIQUE 

We have analyzed 818 pb -1 of e + e~ collisions produced with the Cornell Electron Storage 
Ring at the center-of-mass energy of the f/>(3770) resonance, resulting in about three million 
pairs produced in the CLEO-c detector. We consider candidates where the D° is 
reconstructed as Kg7i°7i°, and the D° is reconstructed in the following four tag modes: D° — >■ 
K + tt~ , D° — > K + 7t~7t°, D° — > K + 7t~7t + 7t~ , and D° — > K + 7t~7t + 7t~7t° (charge congugation is 
implied throughout this paper unless explicitly stated). Reconstructing any of these modes, 
whether it be signal or tag, involves placing requirements on the beam-constrained mass 
m(D) and the energy difference AE, defined as 

m(D) = ^ElUc'-vll? (1) 

AE = E D — -Ebeam (2) 

where -Ebeam is one-half of the center-of-mass energy of the colliding beams, and Ed (pd) is 
the total reconstructed energy (momentum) of the candidate. Both of these quantities rely 
on the fact that Ed would be the same as -Ebeam if the event were perfectly reconstructed, 
and are required to be within 2.5 (3.5) standard deviations of their nominal values for m(D) 
(AE). Figure [1] shows the AE distributions for tagged D° — > Kg7i°7i° candidates that pass 
all other selection criteria. 

We reconstruct neutral pions using photons with E y > 30 MeV and m 77 within 3 standard 
deviations (15 MeV/c 2 ) of the nominal ir° mass 0. We reconstruct Kg — > 7r + 7r~ with m n + n - 
within 2.5 standard deviations (6 MeV/c 2 ) of the nominal Kg mass [7J and with the tc + ty~ 
vertex displaced from the interaction point by at least twice the error on their separation. 
When there are multiple D° — > K^ir°n candidates in a single event we select the one with 
the highest joint probability that both neutral pions are correctly identified. 

After all selection requirements are imposed on both signal and tag-side candidates we 
obtain a combined sample of 1,259 tagged D° — > K^^n events. A Dalitz plot of these 
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FIG. 1: Signal-side AE distribution. The solid line shows the results of fitting with a double 
Gaussian to represent the signal plus a constant background (horizontal dashed line). The vertical 
dashed lines indicate the signal region. 

candidates is shown in Fig. [2j Bose symmetry requires that the decay dynamics remain 
invariant when the two 7r°'s are swapped, hence the Dalitz plot contains two entries per 
event, one for each possible Kgir combination. The horizontal band evident at around 
m 2 (K° s 7i ) = 0.8 GeV 2 /c 4 is due to D° -> K*(892)vr° (if* (892) -> K° s n°), and the most 
striking vertical feature is an empty band at around m 2 (7r°7r°) = 1 GeV 2 /c 4 where the KK 
threshold opens, indicating the presence of destructive interference between the /o(980) and 
some broad underlying im S-wave structure. A study of events both inside and outside the 
signal region indicated in Fig. [1] shows that (7.5 ±0.9)% of the events in this Dalitz plot are 
background. 

We determine the efficiency of our analysis as a function of Dalitz plot coordinate by 
using a GEANT-based Monte Carlo package to generate 1,000,000 simulated D°D° pairs 
where one D is forced to decay to a flavor-tagging mode in proportion to its branching 
fraction while the other D decays to Kg7i°iT uniformly across the Dalitz plot. We pass 
these events through the same analysis code and selection requirements as the data, and 
fit the Dalitz plot of the remaining 63,391 events to a cubic polynomial which has been 
explicitly symmetrized in the two m 2 (KgTT°) variables. We find that the efficiency is well 
modeled by this simple function and is uniform across the Dalitz plot. 

To study backgrounds we use a Monte Carlo-generated sample of D°D° events corre- 
sponding to twenty times the actual integrated luminosity in which the -D's decay to all 
experimentally measured final states with appropriate branching fractions, removing all 
events that contain an actual D° — > KgTT°TT° decay and subjecting the remaining events to 
the same analysis requirements as the real data. The remaining background events vary 
smoothly across the Dalitz plot, preferring the corners of phase space where a 7r° is pro- 
duced at rest and is more easily faked. We find the background shape is well modeled by 
a symmetrized third-order polynomial combined with a small Breit-Wigner component to 
account for u decays from the background mode D° — » Kgu(— » 7r°7) where the u decay 
combines with an unrelated calorimeter shower to yield a second 7r° candidate. 
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FIG. 2: D —> Kgir tt Dalitz plot showing two entries per candidate, one for each possible Kgir 
combination. 

III. DALITZ PLOT ANALYSIS 

The fitting method used to explore the efficiency and background shapes, as well as 
the structure of the D° — > Kgir ^ signal, uses the same unbinned maximum likelihood 
technique described in Ref. js| which minimizes the sum over N events: 



where x and y are the two m 2 (Kgir ) Dalitz variables and V(x, y) is the probability density 
function (p.d.f.) which depends on the event sample being fit: 



The signal p.d.f. is proportional to the matrix element squared, \Ai(x, y)\ 2 , corrected by the 
measured efficiency e(x,y) and the signal fraction, / sig = (0.924 ± 0.009), which is simply 
the complement of the background fraction discussed above, and is fixed in the fit. The 
efficiency, signal, and background contributions are normalized separately: 



N 



£ = -2^1ogP(x n ,y n ) 



(3) 



71=1 




(4) 




(7) 



(6) 



(5) 



providing the overall normalization J V(x, y)dxdy = 1. 
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TABLE I: Details of the resonances included in the Dalitz plot fit. 



Resonance 


Model 


Mass (MeV/c 2 ) 


Width (MeV/c 2 ) 


7T7T pole 


Complex Pole 


470 - i220 




K*(892) 


Breit-Wigner 


896 


50.3 


/ (980) 


Flatte a 


965 




/ 2 (1270) 


Breit-Wigner 


1275.1 


185.0 


/o(1370) 


Breit-Wigner 


1350 


265 


if|( 1430) 


Breit-Wigner 


1432.4 


109 


/o(1500) 


Breit-Wigner 


1505 


109 


if* (1680) 


Breit-Wigner 


1717 


322 



a g„„ = 406 MeV/c 2 , g KK /9™ = 2 

Several models for fitting the data have been explored, and three of these are described 
here. All include the if* (892), K%(1430), and if* (1680) as intermediate if^vr resonances 
and an / 2 (1270) intermediate 7r°7r° resonance. To model the S-wave contribution to 7r°7r° we 
add an / (980), described with a Flatte parameterization |9j, to a ttti complex pole around 
the mass of the <r(500), and we include one of the following three variations: /o(1370) 
(Model 1), / (1500) (Model 2), and both / (1370) and / (1500) (Model 3). Each includes 
a non-interfering Kg — > tt°tt° contribution to account for D° — > K S K S . A fourth model, 
in which the broadest S-wave features were replaced by a simple non- resonant component, 
was rejected since the fit fraction of the non-resonant component exceeded 80%. Table d 
summarizes the parameters used for the intermediate resonances and the functions used to 
represent them. 

Table [III summarizes the Dalitz plot fit results for the three models described above. In 
all cases we choose the amplitude (phase) of the if* (892) to be 1 (0) respectively, which 
effectively defines the complex coordinate system for the other resonances included in the 
fit. The x 2 values shown are calculated post-fit by dividing the Dalitz plot into discrete bins 
and comparing the data with the average p.d.f. in each bin. The fit fraction (FF) for each 
component is calculated by integrating its contribution across the Dalitz plot and dividing 
by the integral of the coherent sum of all components: 



FF l = 



Ivv 


Oil 


2 


Ivv 




a J 


e icf>j 


2 



(8) 



where aj ((f) j) are the amplitude (phase) of the j-th component. The total fit fraction for a 
given model is the sum of the fit fractions for the individual components of that model, and 
a total fit fraction greater than unity indicates the presence of destructive interference on 
the Dalitz plot. 

While the quality of the fit projections of the three models are comparable, we choose 
Model 1 as our preferred fit since it has both a reasonable \ 2 (unlike Model 2) and a total 
fit fraction close to unity (unlike Model 3). Figures [3] and H] show the projections of both 
data and p.d.f. for Model 1. In both figures the data are represented by the points with 
error bars, the total fit by the thick solid line, the background contribution by the thin solid 
line, the tttt S-wave contribution by the dashed line, and the if* (892) contribution by the 
broken solid line. 

The results in Table [III show a large model dependence of the fit fractions and phases of 
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FIG. 3: m 2 (Kir) projections for Model 1. The 
data are represented by the points with error 
bars and the total fit by the thick solid line. 
The other fit components are discussed in the 
text. 



FIG. 4: m 2 (7T7r) projection for Model 1. The 
data are represented by the points with error 
bars and the total fit by the thick solid line. 
The other fit components are discussed in the 
text. 



both the /o(980) and the tttt pole while the overall fit projections remain largely unchanged. 
We interpret this as an indication that, while the overall shape and phase evolution of the 
7T7T S-wave component of our fits are reasonable, the details of the underlying physics model 
used to describe this S-wave component are not well determined. To test the hypothesis that 
the combined structure of the S-wave components is meaningful while the specific details 
of these models are not, we combine all of the individual spin-zero 7r°7r° components into a 
single tttt "S-wave object" and calculate the fit fraction of this object as a whole: 



Sxrp 




<t>k 


2 


Ivv 


EoH a 3 ei4>i 


2 



(9) 



The results summarized in Table III II show that the tttt S-wave fit fraction is now consistent 
among models, as are the total fit fractions which are now near 100%. 

The fact that the underlying physics is not well represented by a simple combination of 
resonances is not surprising given the results of other analysis efforts where "a" like objects 
have turned out to be the result of dynamical effects that require a more sophisticated 
modeling approach |Tlj ]. In spite of this shortcoming, our simple model does capture the 
striking overall feature that there is significant destructive interference at m 2 (7r°7r°) around 
1 GeV 2 /c 4 , and that the overall ir ir S-wave component accounts for about a quarter of all 
D° — > KgTT°iT decays. These results are consistent with the studies of D° — > K^n+ir' from 
BaBar and Belle in which the total 7T7t S-wave plus non-resonant fit fractions were found to 
be 27.2% and 26.1% respectively 

To test the stability of our results as selection requirements are varied, we tightened 
(loosened) these requirements by half a standard deviation resulting in 1,099 (1,361) candi- 
dates in the Dalitz plot. In each case the signal fraction / sig , efficiency shape e(x,y), and 
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TABLE II: Fit fraction (%), amplitude, and phase (°) fit results for Models 1 through 3. The errors 
are statistical only. 







ivioaei i 


ivioaei z 


ivioaei 


X 2 /N dof 


1 Q Q /7 


zLzL 7 /7 


1 4 Q 


Total FF 


IZZzto 


1 9D-I-1 1 


9c;9_i_QQ 
ZOZztoo 




FF 


2.7±1.4 


3.7±1.8 


4.5±2.2 


7T7T pole 


a 


0.67±0.16 


0.91±0.20 


0.99±0.23 




4> 


140±17 


119±22 


39±17 




FF 


10.5±2.1 


12.1±2.4 


18.4±4.3 


/o(980) 


a 


l.7l±0.l7 


2.13±0.20 


2.59±0.24 




<P 


35.2±9.9 


65±11 


44.8±7.9 




FF 


25.7±5.l 


(n/a) 


81±24 


/o(1370) 


a 


5. f 2±U.5o 


(n/a) 


11.0±1.0 




<f> 


o4U.o±0.t> 


(n/aj 


ID.ozho.o 




FF 


(n/aj 


99 q_i_« n 




/o(1500) 


a 


I'm /al 


1 1 7_i_i c 


9fl Q-I-/L n 




4> 


(n/ aj 


1D±1Z 


901 /14-c n 

Zol .4±o.U 




FF 


z.4o±u.yi 


1 9 n_i_ Q Q 


O.OztZ.D 


/ 2 (1270) 


a 


l.O I ICU.ZO 




9 o«-i-n ^ 






9K9-I-1 a 
ZoZztlo 


9 94-fi ^ 


Q/in Q4-8 Q 

o4u.y±o.y 




FF 


65.6±5.3 


48.6±5.9 


50.2±7.8 




a 


I (nxeaj 


l (nxeaj 


l (nxeaj 




<f> 


u ^nxeci j 


r» /Xj vor |\ 

u ^nxeQ j 


u ^nxeo ) 




FF 


0.49±0.45 


1.9±1.2 


1.45±0.82 


#3(1430) 


a 


0.43±0.18 


0.98±0.29 


0.85±0.24 




4> 


141±28 


191±16 


159±15 




FF 


11.2±2.7 


15.2±4.3 


13.5±3.5 


K*(1680) 


a 


5.65±0.77 


7.6±1.2 


7.07±0.82 




<f> 


55±11 


45±12 


18.7±8.9 




FF 


3.46±0.92 


3.30±0.96 


2.48±0.95 




a 


0.281±0.037 0.318±0.044 0.272±0.050 






contributes incoherently 



TABLE III: Dalitz plot fit fractions (%) with a coherent tttt S-wave object. 



Model 1 Model 2 Model 3 



7T7T S-wave 


28.9±6.3 


19.5±5.7 


30±12 


K*(892) 


65.6±5.3 


48.6±5.9 


50.2±7.8 


Kl (1430) 


0.49±0.45 


1.9±1.2 


1.45±0.82 


#*(1680) 


11.2±2.7 


15.2±4.3 


13.5±3.5 


/ 2 (1270) 


2.48±0.91 


12.9±3.3 


6.8±2.5 




3.46±0.92 3.30±0.96 2.48±0.95 



Total 112.1±8.8 101±10 105±15 
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TABLE IV: Dalitz plot fit results for Model 1. The first error is statistical and the second represents 
variations on analysis requirements. 

Component Fit Fraction (%) 

7T7T S-wave 28.9 ± 6.3 ± 3.1 

K*(892) 65.6 ± 5.3 ± 2.5 

K 2 *(1430) 0.49 ± 0.45 ± 0.23 

K*(1680) 11.2 ± 2.7 ± 2.5 

/ 2 (1270) 2.48 ± 0.91 ± 0.78 

K Q S 3.46 ± 0.92 ± 0.66 



background shape B(x, y) were remeasured and fixed in the fit to the Dalitz plot, which used 
Model 1 as described above. Sensitivity to uncertainties in identifying and reconstructing 
7T°s was studied by tightening the nominal pull mass requirement to ±2.5 standard devi- 
ations, resulting in 1,150 events in the Dalitz plot, and by eliminating all ir°s having one 
or more photons in the endcap calorimeter (cos^ > 0.81), resulting in 1,023 events in the 
Dalitz plot. In the same way, sensitivity to the signal fraction was studied by varying / s ; g 
by ±3 standard deviations from its nominal value; sensitivity to the shape of the efficiency 
function was studied by replacing the fitted efficiency shape parameters with those corre- 
sponding to a uniform efficiency e(x,y) = 1; and sensitivity to uncertainty in the shape of 
the background was studied by replacing the nominal background shape, determined us- 
ing Monte Carlo-generated data, with a background shape determined by analyzing single 
tagged D° — >■ Kg7T°7i° data selected from a sideband region in the AE vs. m(D) plane. 

In all cases described above, the parameters derived from fitting the resulting Dalitz plot 
with Model 1 were within 1 standard deviation of the parameters found in the preferred 
fit. For each parameter we choose the maximum variation from the preferred fit value to 
represent the systematic uncertainty in that parameter. The resulting systematic errors for 
the fit fractions from Dalitz plot Model 1 are given in HVl 



IV. BRANCHING FRACTION ANALYSIS 

The extraction of the D° — > Kgii 7T branching fraction is a straightforward extension of 
the Dalitz analysis described above; we need only to quantify the yields and efficiencies for 
observing the signal D° — > Kg7i 7r° events as a function of the tag mode. More specifically, 
to extract the D° — > K^tt ?! branching fraction using double-tagged events where one D 
decays to Kg7i°7i° and the other D decays to tag mode r, we need to evaluate 

B{D K s ii 7T ) T = — — — , (10) 

61 v £,0750 £V 

where e T DT is the efficiency for the detector and analysis to select an event where the D° 
decays to i^7r°7r° and the D° decays to tag mode r, iV£, T is the measured yield of these 
events, N D oj=j is the number of D°D° pairs produced, B T is the branching fraction of the 
tag mode, a T is a small correction to account for the quantum correlation between the two 
sides of the event, and the factor of 2 accounts for the fact that we include charge-conjugate 
decays. 
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TABLE V: The yields, efficiency, and quantum correction for each tag mode. 



Tag Mode N T DT e T DT (%) 2N D0 - m B T a T 

K + ir- 247±17 9.29±0.41 229050±600 1.130±0.041 
K+tt-tt 500±25 5.30±0.26 809700±1700 1.099±0.027 
K+7T-TT+TT- 358±28 6.04±0.30 449500±1100 1.080±0.028 



The yield of D° — > KgTT°TT° signal events for each tag mode, N^ T , is determined by fitting 
the AE distribution for the tagged events to a double Gaussian signal plus a flat background 
as illustrated in Fig. [TJ The yields for all tag modes are summarized in Table IVl 

The efficiency for detecting a D° decaying to K^^n when the D° decayed to tag mode 
t was determined by generating a large sample of simulated D°D° events for each case 
and measuring the fraction of these events that pass the event selection requirements. The 
numbers are summarized in Table N] and include an additional factor of e c °o r = 0.94 ± 0.02 
for each 7r° to correct for a known systematic difference between data and Monte Carlo. 

The factor 2N DO j^oB T in the denominator of Eq. [10] is number of the P's that decayed to 
tag mode r in the CLEO-c detector. This is determined from the efficiency-corrected single- 
tag yields obtained by a recently published CLEO-c a naly sis that used the same data sample 
and tag modes to study semileptonic charm decays [12j. The numbers are summarized in 
TableED 

The correction factor a T is needed because the signal-side D decays to a CP-even eigen- 
state and the tag-side D decays to a flavor eigenstate of mixed CP, allowing interference 
between favored and suppressed decay paths to modify the observed decay rate at the five 
percent level. The formalism is discussed in Refs. |3[ and [HI where a T is defined as 

a T = l + 2rcos5 T + R T ws + y = l + R T ws + A T QC (11) 

and S T is the strong phase difference between the amplitudes of D° and D° decaying to tag 
mode r, R^s * s ^ ne relative rate of D° and D° decaying to tag mode r, and y is the standard 
D°D° mixing parameter while the values of and Aq C for the D° — > K + ir~ tag mode 
are obtained from the Particle Data Group [3], and the values for the D° — > K+n^n and 



D° — > K + 7i 7r + 7r tag modes are obtained from Ref. 13j|. Since PJy 5 and Aq C are not 
known for D° — > i^ + 7r~7r + 7r~7r°, this tag mode is not used to extract the D° — > Kg7r°7T 
branching fraction. The overall correction factor a T for each tag mode is reported in Table PVT. 
We calculate the D° K^n n° branching fraction for each tag mode and show the results 
in Table [VB 

To study the systematic uncertainty associated with these measurements we varied the 
way both the signal and background are parameterized when extracting signal yields, in 
all cases performing the same procedure on both data and Monte Carlo. Other analysis 
variations were also explored, including the tightening and loosening of reconstruction re- 
quirements as described in the previous section. These variations were combined in quadra- 
ture with the uncertainties in a, e c °S r ■, Nst/^st, tracking, Kg efficiency, and the statistical 
uncertainty in measuring Edt, to give the total systematic uncertainty for each tag mode. 
Table IVII lists the results from our branching ratio calculations and gives our final result by 
performing a weighted average across the three modes. 
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TABLE VI: Summary of B(D° — > K^^ir ) measurements. 

Tag Mode B(D° K%ir°ir ) (%) 
K + ir° 1.030 ± 0.069 ± 0.094 
K + ir~ir 1.061 ±0.054 ±0.110 
K+tt-tt+tt- 1.099 ±0.084 ±0.115 
Average 1.059 ± 0.038 ± 0.061 



V. CONCLUSION 

We have performed a Dalitz plot analysis of the decay mode D° — > Kg7i°7v° using the 
full CLEO-c data set and have measured the total branching fraction of this mode to be 
(1.059 ± 0.038 ± 0.061)%. We find that the combined 7r°7r° S-wave contribution to our 
preferred fit is (28.9±6.3±3.1)% of the total decay rate while D° ->• F*(892)°7r° contributes 
(65.6 ± 5.3 ± 2.5)%. 
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